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ABSTRACT 

We describe a measurement of the angular power spectrum of anisotropics in the Cosmic Microwave 
Background (CMB) from 0.3° to 10° from the North American test flight of the Boomerang ex- 
periment. Boomerang is a balloon-borne telescope with a bolometric receiver designed to map CMB 
anisotropies on a Long Duration Balloon flight. During a 6-hour test flight of a prototype system in 1997, 
we mapped > 200 square degrees at high galactic latitudes in two bands centered at 90 and 150 GHz with 
a resolution of 26 and 16.6 arcmin FWHM respectively. Analysis of the maps gives a power spectrum 
with a peak at angular scales of ^ 1 degree with an amplitude ~ lOfiKcAiB- 

Subject headings: cosmology: Cosmic Microwave Background, anisotropy, measurements, power 
spectrum 



1. INTRODUCTION 

Measurements of Cosmic Microwave Background 
(CMB) anisotropies have the potential to reveal many 
of the fundamental properties o f the universe (e.g. see 



Kamionkowsi fc Kosowsky 1999 and references therein) 



Since the mea surement of the larg e scale anisotropies by 
COBE-DMR ([Wright et al. 1994|, [Bennett et al. 1996| , 



etc.), many ground-based and balloon-borne experiments 
have continued to develop the technology necessary to pro- 
duce accurate measurements of CMB structure on smaller 
angular scales. Recent results from these experiments have 
shown evidence for the existence of a peak in the power 
spectrum of CMB fluctu a tions at a multipole, £ ^ 200 
([Ncttcrfield et . al. 19971 [dc Oliveira-Costa et al. 1998 



Wilson et al. 1999|, iTorbet ct al. 19991 , [Coble et al 



1999, Miller ct al. 1999). In this paper, we present re- 



sults from the test flight of Boomerang that constrain 
the position and amplitude of this peak. In fact, the re- 
sults from this data set alone can be used to constrain cos- 
mological models and provide evidence for a flat universe, 
Q — ilM + — 1 (see the companion paper, Melchiorri, 



et al. 1999). 

Boomerang is a millimeter- wave telescope and receiver 
system designed for a Long Duration Balloon (LDB) flight 
from Antarctica of up to two weeks duration. The data 
presented here are from a prerequisite test flight of the 
payload in August, 1997, which lasted 6 hours at float alti- 
tude. A prototype of the LDB focal plane was used to test 
the experimental strategy of making a map and measur- 
ing the angular power spectrum of CMB anisotropy with 
a slowly-scanned telescope and total-power bolometric re- 



2. INSTRUMENTATION 

Boomerang employs bolometric detectors cooled to 
0.3 K, coupled to an off-axis telescope and operated as 
total power radiometers. Electronic modulation and syn- 
chronous demodulation of the bolometer signals provide 
low frequency stability of the detector system. The opti- 
cal signal is modulated by slowly scanning the telescope 
(~ 1 deg/s) in azimuth. A descrip tion of the instrument 
can be found in Masi ct al. 1999a. 
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The Boomerang telescope consists of three off-axis, 
aluminum mirrors: an ambient temperature f .3 m primary 
mirror and two smaller (15 cm) 2 K reimaging mirrors. 
The telescope is shielded from ground radiation by exten- 
sive low emissivity baffles. Incoming radiation is reflected 
by the primary mirror into the entrance window of the 
cryostat and reimaged by the 2 K optics onto the 0.3 K fo- 
cal plane. The tertiary mirror is positioned at an image of 
the primary mirror and forms a throughput-limiting Lyot 
stop illuminating the central 85 cm of the primary mirror. 
The inside of the 2 K optics b ox has a thick coating of 
absorbing material (Bock 1994) to intercept stray light. A 
calibration lamp mounted in the Lyot stop of the optical 
system is pulsed for ^ 0.5 s every 15 minutes to produce 
a stable, high signal-to-noise ratio ( ^ 1000) calibration 
transfer. 

For the test flight, the foc al plane contained silicon ni- 
tride microm esh bolometers (Mauskopf et al. 1997, Bock 



et al. 1998) fed by multi-mode feed horns with channels 
at 90 and 150 GHz. Fourier transform spectroscopy of 
each channel determined the average band centers and ef- 
fective band widths for sources with a spectrum of CMB 
aniso tropy ( Table 1). A long d uration ^He cryostat (Mas] 
et al. 1998| , [Masi et al. 1999b ) maintained the focal plane 
at 0.285 ± 0.005 K during the flight. The bolometers were 
biased with an AC current and the signals were demodu- 
lated, amplified and passed through a 4-pole Butterworth 
low-pass filter with a cut-off frequency of 20 Hz for the 
150 GHz channels and 10 Hz for the 90 GHz channels. 
The bolometer DC levels were removed with a single-pole, 
high-pass filter with a cut-on frequency of 16 mHz and re- 
amplified before digital sampling at 62.5 Hz with 16 bit 
resolution. 

The attitude control system included a low-noise two- 
axis flux-gate magnetometer, a CCD star camera and three 
orthogonal (azimuth, pitch, roll) rate gyroscopes. The az- 
imuth gyroscope drove the azimuth feedback loop, which 
actuated two torque motors, one of which drove a flywheel 
while the other transferred angular momentum to the bal- 
loon flight line. The gyroscope was referenced every scan 
to the magnetometer azimuth to remove long time scale 
drifts. Pendulations of the pay load were reduced by a oil- 
flUed damper mounted near the main bearing and sensed 
by the pitch and roll gyroscopes outside of the control loop. 
The CCD frame had 480 x 512 pixels, with a resolution 
of 0.81 arcmin/pixel (azimuth) and 0.65 arcmin/pixel (el- 
evation) . The limiting magnitude of stars used for point- 
ing reconstruction was my ~ 5. The coordinates of the 
two brightest stars in each CCD frame were identifled by 
on-board software and recorded at 5 Hz, for post-flight 
attitude reconstruction. 

The primary mirror and the receiver were mounted on 
an inner frame which could be rotated to position the ele- 
vation boresight of the telescope between 35° — 55°. During 
observations, the inner frame was held flxed while the en- 
tire gondola slowly scanned in azimuth. The scan speed 
was set to the maximum allowed by the beam size and 
bolometer time constants. The window functions for the 
90 and 150 GHz channels in this flight were rolled off at 
i ^ 300 and 500, respectively, by the 26 and 16.5 arcmin 
beams and were limited to £ 20 by the system 1 /f noise. 

3. OBSERVATIONS 



The instrument was launched at 00:25 GMT on August 
30, 1997 from the National Scientific Balloon Facility in 
Palestine, Texas. Observations started at 3:50 GMT and 
continued until the flight was terminated at 9:50 GMT. 
Sunset was at 01:00 GMT and the moon rose at 10:42 
GMT. The average altitude at float was 38.5 km. We ob- 
served in three modes: ~ 4.5 hours of CMB scans, 40 min- 
utes of calibration scans on Jupiter, and four ten-minute 
periods of full-sky rotations of the gondola. The CMB 
observations consisted of smoothed triangle wave scans in 
azimuth with a peak-to-peak amplitude of 40 degrees cen- 
tered on due South at an elevation angle of 45 degrees. 
The scan rate was ^ 2.1 deg/s in azimuth (1.4 deg/s on 
the sky) in the linear portion of the scan (63% of the pe- 
riod). These scans, combined with the earth's rotation, 
covered a wide strip of sky, from —73° < RA < 23° and 
-20° < DEC < -16°. 

4. DATA REDUCTION 

We obtain a first order reconstruction of the attitude 
from the elevation encoder and magnetometer data. Us- 
ing this, the stars recorded by the CCD camera are identi- 
fied. The data from the gyroscopes are used to interpolate 
between star identifications. The offset between the mi- 
crowave telescope and the star camera is measured from 
the observation of Jupiter. The attitude solution obtained 
in this way is accurate to ~ 1 arcmin rms. 

The results presented in this paper are primarily from 
the highest sensitivity 150 GHz channel. Data from a 
90 GHz channel are presented as a systematic check. The 
data at float consist of '--^ 1.4 million samples for each 
bolometer. These data are searched for large (> 5(t) de- 
viations such as cosmic ray events and radio frequency 
interference and flagged accordingly. Flags are also set 
for events in the auxiliary data such as calibration lamp 
signals and elevation changes, as well as for the different 
scan modes. Unflagged data, used for the CMB analysis, 
are 54% of the total at 150 GHz and 60% of the total at 
90 GHz. We deconvolve the transfer functions of the read- 
out electronics and the bolometer thermal response from 
the data, and apply a flat phase numeric filter to reduce 
high frequency noise and slow drifts. The time constants 
of the bolometers given in Table 1 are measured in flight 
from the response to cosmic rays and to fast (18 deg/s) 
scans over Jupiter. 

Table 1. BOOMERANG 1997 instrument parameters 
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The primary calibration is obtained from scans of the 
telescope over Jupiter, which were performed from 3:32 
U.T. to 3:53 U.T. (1.2 deg/s azimuth scans). Jupiter was 
also re-observed later during regular CMB scans (from 5:58 
to 6:18 U.T.) and during full sky rotations. Jupiter was 
rising during the first observation (from an elevation of 
36.9° to 38.2°) and setting during the second one (from 
an elevation of 40.1° to 39.0°). The deconvolved and fil- 
tered data from these scans are triangle-interpolated on a 
regular grid centered on the optical position of Jupiter to 
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make a beam map. The beams are symmetric with minor 
and major axes equivalent within 5%. Sohd angles given 
in Table 1 are computed by integration of the interpolated 
data. The beam map is used to derive the window func- 
tion, Wi, and an overall normalization. 

We use a brightness temperature of T^ff ~ 17 3 K for 
both the 90 and 150 GHz bands ( |Uhch 198l| , |GrifRn| 
1986, iHildcbrand 1985 ) and assume an uncertainty of 5% 
( Goldin et al. 1997 ). Wc also assign a 5% error in the con- 
version from brightness temperature to CMB temperature 
due to statistical noise in the measurements of the band 
pass. Errors in the determination of the receiver transfer 
function are largest at the highest temporal frequencies 
where signals are attenuated by the bolometer time con- 
stants. These errors affect the calibration of the window 
function and the CMB power spectrum measurements at 
the largest multipoles. These should cancel because the 
beam maps of Jupiter are made from scans at the same 
scan speed as the CMB scans. In simulations, a 15% error 
in the determination of the bolometer time constant pro- 
duces a maximum error < 1.5% in the normalization and 
^ 6% in the value of at £ < 300 near scan turnarounds 
where the scan speed was only 1 deg/s. The final precision 
of the calibration including errors in the temperature of 
Jupiter, Tcff, the transfer function used for deconvolution, 
the measurement of the filter pass bands and the beam 
solid angle is 8.1% and 8.5% at 150 GHz and 90 GHz, 
respectively. 

Signals from the internal calibration lamp are used to 
correct for a slow linear drift in detector responsivity dur- 
ing the flight, mainly due to variation (0.29-0.28 K) in the 
base temperature of the fridge. The change in response 
to the calibration lamp from the beginning to the end of 
the flight is -1-9%, -1% at 90 GHz, and -1-8%, -1.5% at 
150 GHz. 

5. MAP MAKING AND POWER SPECTRUM ESTIMATES 

Current and future CMB missions require new meth- 
ods of analysis able to incorporate the effects of corre- 
lated instrument noise and new implementations capa- 
ble of processing large data sets (Bond et al. 1999). 
The analysis of the Boomerang test flight data provides 
a test of these methods on a moderate size (~ 25,000 
pixels) data set. The calibrated time stream data are 
processed to produce a pixelized map, and from this a 
measurement of the angular power spectrum, using the 
MADCAP software package of BorriU ( 1999a, 1999b) (see 



http://cfpa.berkeley.edu/^borrill/cmb/madcap.html) on 
the Cray T3E-900 at NERSC and the Cray T3E-1200 at 
CINECA. 

Noise correlation functions are estimated from the time 
stream under the assumption that in this domain the sig- 
nal is small compared to the noise. For the 150 GHz chan- 
nel, the maximum likelihood map with 23,561, 1/3-beam 
sized (6.9') pixels is calculated from the noise correlation 
function, the bolometer signal (excluding flagged data and 
all data within 2 degrees of Jupiter) and the pointing infor- 
mation using the procedures described in Wright (1996), 
Tegmark (1997) and Ferreira & Jaffe (1999). The maxi- 
mum likelihood power spectrum is estimated from the map 
using a Newton-Raphson iterative maximization of the Ci 
likelihood function following Bond, Jaffe & Knox (1998) 
(BJK98) and Tegmark (1998). The power spectrum analy- 



sis requires approximately 12 hours on sixty-four T3E-900 
processors. 

Table 2: Orthogonalized power spectrum measurement 
from boomerang/na150 ghz channel map with with 6 ar- 
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Fig. 1. — Power spectrum of the Boomerang 150 
map with 6 arcminute pixelization. The solid curve 
is a marginally closed model with {i}b,^M ,^A,ns , h) = 
(0.05, 0.26, 0.75, 0.95, 0.7). The dotted curve is standard CDM 
with (0.05,0.95,0.0,1.0,0.65). Dashed curves are open and 
closed models with fixed r 2totai = 0.66 and 1.55, res pectively 
(see the companion paper, A. Melchiorri, et al. 1999). 



Due to finite sky coverage, individual Ci values are 
not independent, and we estimate the power in eight 
bins over 25 < ^ < 1125, binning more finely at low £ 
than high. In each bin we calculate the maximum likeli- 
hood amplitude of a flat power spectrum in that bin, i.e. 
£{£ + l)Ci ^constant (Table 2). Given our insensitivity 
to low £ signal, and in order to ameliorate the effects of 
low-frequency noise (which translates into \ow-£ structure 
for this scanning strategy), we marginalize over power at 
^ < 25 and diagonalize the bin correlation matrix using 
a variant of techniques discussed in BJK98. Bin to bin 
signal correlations are small but non-negligible, with each 
top-hat bin anticorrelated with its nearest neighbors at 
approximately the 10% level. 

We calculate likelihood functions for each bin using the 
offset log-normal distribution model of BJK98. We quote 



errors for the power spectrum from ±1ct errors assuming 
simple Gaussian likelihood functions and from 68% con- 
fidence intervals for the full likelihood functions. Plots 
of the likelihood functions are presented in the compan- 
ion paper (Melchiorri, et al. 1999) and full spectral data 
including information on the shape of the likelihood func- 
tion (Bond, Jaffe & Knox 1998b) and appropriate window 
and filter functions (Kn ox 1999) will be available at t he 
Boomerang web site ( | http://boom. physics. ucsb.edu or 



http: / /oberon.romal .infn.it /boom ) . In Figure 1, we show 
the calculated power spectrum for the 150 GHz channel 6 
arcminute pixel map. 

6. SYSTEMATICS CHECKS AND FOREGROUNDS 

The scan strategy and analysis allow for a variety of 
checks for systematic effects in the data. For all system- 
atic tests we produce small (~ 4000 pixel) maps with 16' 
pixels which can be generated and analyzed quickly (~ 30 
minutes on four T3E-900 processors). Monte-Carlo simu- 
lations show the effects of the 16' pixelization are negligible 
for values of ^ < 400. Compared to the level of statisti- 
cal noise, we find that the derived power spectrum from 
the 150 GHz channel is insensitive to (i) the size of the 
excluded region around Jupiter, (ii) the choice of spectral 
shape within the multipole bins, and (iii) the method of 
marginalizing over low £ signals. 

We also analyze maps made from combinations of data 
that we expect to produce null power spectra in the ab- 
sence of spurious sources of noise power: (i) data from a 
dark detector and (ii) the difference between left and right- 
going scans at 150 GHz. In each case (Table 3), we find a 
power spectrum consistent with zero. 

Table 3. Power Spectra, £{£ + l)Ci/2ir, from BOOMERANG sys- 
tematic tests. Error bars are itlcr, assuming Gaussian likelihood 
functions. Unlike the data in Table 2 and Figure 1, these spectra 
are calculated with 16 arcminute pixels and have not been orthog- 
onalized. The (L-R) difference analysis in column 5 has not been 
corrected for the effects of correlated noise, which leads to negative 
power estimates in the high £ bins. 
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The sky strip observed is extended in Galactic latitude 
from 6 ~ 15° to 5 ^ 80°. From IRAS/DIRBE map ex- 
trapolation (Schlegel et al. 1998) dust emission is not ex- 
pected to produce significant contamination in this region 
of the sky. We find an amplitude for a flat power spectrum 
from 25 < e < 400 of £{£ + l)Cf/27r = (3580 ± 706)AiK2 
and (2982 ± 702)^K^ from the half of the map near the 



galactic plane (15° < b < 45°) and the higher-latitude half 
(45° < 5 < 80°), respectively. Because of sky rotation, this 
also corresponds to separately analyzing the data from the 
first half and second half of the flight. 

Jupiter is more than 1000 times brighter than degree- 
scale anisotropy in the CMB and is in the middle of the 
lowest latitude half of the map. We test for sidelobe con- 
tamination by varying the size of the circular avoidance 
zone around Jupiter from (j) = 0.5 — 3 degrees for the CIVIB 
maps, and find no significant change in the power spectrum 
as long as we remove data within ^ 1 degree. 

7. DISCUSSION AND CONCLUSIONS 

Data obtained in ~ 4.5 hours of CMB scans during the 
Boomerang test flight have been analyzed to produce 
maps of the millimeter- wave sky at 90 and 150 GHz, with 
26 and 16 arcmin FWHM resolution, respectively. The 
instrument was calibrated using observations of Jupiter. 
The experiment employed a new scan strategy and detec- 
tor technology designed to give maximum coverage of an- 
gular scales. The data analysis implements new techniques 
for making maximum likelihood maps from low signal-to- 
noise time stream data over large numbers of pixels. The 
power spectrum Ce obtained from the maps extends from 
^ ~ 50 to ^ - 800 and shows a peak at £ - 200. 

Since the BoOMERANG/NAtest flight in August, 1997, 
we have obtained > 200 hours of data from the LDB 
flight of Boomerang (Boomerang /LDB), carried out 
in Antarctica at the end of 1998. The focal plane for the 
LDB flight contained 2, 6, and 3 detectors at 90, 150, and 
240 GHz, each with better sensitivity to CMB, faster time 
constants, and lower 1 /f noise than the best channel in the 
BooMERANG/NAflight, as well as 3 detectors at 400 GHz 
to provide a monitor of interstellar dust and atmospheric 
emission. The results from this flight will be reported else- 
where. 
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